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We study the effect of a strong static homogeneous electric field on the highly excited rovibrational 
levels of the LiCs dimer in its electronic ground state. Our full rovibrational investigation of the 
O ' system includes the interaction with the field due to the permanent electric dipole moment and 

D I the polarizability of the molecule. We explore the evolution of the states next to the dissociation 

, threshold as the field strength is increased. The rotational and vibrational dynamics are influenced 

by the field; effects such as orientation, angular motion hybridization and squeezing of the vibrational 
motion are demonstrated and analyzed. The field also induces avoided crossings causing a strong 
mixing of the electrically dressed rovibrational states. Importantly, we show how some of these 
highly excited levels can be shifted to the continuum as the field strength is increased, and reversely 
' how two atoms in the continuum can be brought into a bound state by lowering the electric field 

^ I strength. 

PACS numbers: 32.60.+i,33.20.Vq,33.80.Gj 
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O ■ I. INTRODUCTION 
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The recent availability of cold to ultracold polar dimers in the vibrational and rotational ground state of their 
singlet electronic ground potential P, 0, represents a breakthrough towards the control of all molecular degrees 
of freedom, i.e., the center of mass, electronic, rotational and vibrational motions, and towards the ultimate goal of 
obtaining a polar condensate. The experimental achievements have been accompanied by significant theoretical efforts 
to understand the intriguing physical phenomena expected for ultracold polar quantum gases due to their anisotropic 
and long-range dipole-dipole interaction. In particular it has been analyzed how external fields can control and 
manipulate the scattering properties^ JBL Q , and the chemical reactions dynamics Q , or how to use them as tools 
. for quantum coniputational devices [3, [3, • Different approaches to achieve cold and ultracold molecules have been 
Tj" ' explored [ill [13) [3, [3 [IB| • For further work and references we refer the reader to the comprehensive reviews [H, [13] ■ 
, The most widespread techniques to produce ultracold polar molecules are the photoassociation of two ultracold 
' atoms [3, [3] , and the tuning of the atomic interactions via magnetically induced Feshbach resonances [1^ . Al- 
\ ternative pathways explore the ability to manipulate the interaction between atoms by inducing optical Feshbach 
. resonances. Based on the same principle as the magnetically induced Feshbach resonances, they appear when two 
^ ■ colliding ultracold atoms are coupled to a bound state of the corresponding molecular system by using a radiation 
• ^ . field. Initially, there were several theoretical proposals to obtain these resonances with the help of radio-frequency, 
' static electric, and electromagnetic fields [2l|, [2^, [2^ [2^ . In addition, it has been demonstrated that a combination of 
^ ', a magnetic and static electric field can induce Feshbach resonances in a binary mixture of Li and Cs atoms [25, 2(|, 
■ - - I and that a suitable combination of these two fields can tune the relevant interaction parameters, such as the width and 



open-channel scattering length in these resonances [27[. The existence of this optically induced resonances has been 
experimentally proved for different atomic species by tuning the laser frequency near a photoassociation resonance 

[11, m [13, ill. 

Within the above experimental techniques, the molecules are usually in a highly excited vibrational level close to the 
dissociation threshold of an electronic state. These vibrational states are exposed to external fields and most of their 
overall probability is located at the outermost hump of their probability density. In the present work, we investigate 
the few last most weakly bound states of the X^S]+ electronic ground state of a polar molecule in a strong static electric 
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field. We perform a full rovibrational investigation of the field-dressed nuclear dynamics including the interaction of 
the field with the molecular electric dipole moment and polarizability. The LiCs dimer is a prototype system and will 
be used here. This choice is based on the experimental interest in this system and the availability of its molecular 
polarizabilities . It completes our previous investigations on the effects of an electric field on this system, where we 
have analyzed the rovibrational spectrum, the radiative decay properties, and the formation of these ultracold dimers 
via single-photon photoassociation from the continuum into its electronic ground state [13, [3, [13, [sH] • Specifically, 
we analyze the binding energies and the expectation values, (cos 6'), (J^) and (i?), of states lying in the spectral region 
with binding energies smaller than 0.28 cm~^ and vanishing azimuthal quantum number in the very strong field 
regime. We demonstrate that both the rotational and vibrational dynamics are significantly affected by the field. 
Indeed, the vibrational motion is squeezed or stretched to minimize the energy, depending on the rotational degree 
of excitation and the field strength. At such strong fields the nuclear spectrum exhibits several avoided crossings 
between energetically adjacent states, which lead to a strongly distorted rovibrational dynamics. The latter might 
be directly observable when imaged by photodissociation experiments. Beyond this, (magnetically induced) avoided 
crossings in Cs2 have been used to construct a molecular Stiickelberg interferometer [33 . In addition, we show that by 
tuning the electric field strength a dissociation channel is opened, i.e., a weakly bound molecular state with low- field 
seeking character is shifted to the atomic continuum by increasing the field strength. Of course, the reverse process 
is also possible, and two free atoms can be brought into a molecular bound state by lowering the field strength. This 
might be of interest to control the coUisional dynamics of the atomic/molecular cold gas by using either very strong 
static (micro-) electric fields or strong quasistatic, i.e., time-dependent fields. 



II. THE ROVIBRATIONAL HAMILTONIAN 



We consider a heteronuclear diatomic molecule in its electronic ground state exposed to a homogeneous and 
static electric field. Our study is restricted to a non-relativistic treatment and addresses exclusively a spin singlet 
electronic ground state, and therefore relativistic corrections can be neglected. We assume that for the considered 
regime of field strengths perturbation theory holds for the description of the interaction of the field with the electronic 
structure, whereas a nonperturbative treatment is indispensable for the corresponding nuclear dynamics. In addition, 
we take into account the interaction of the field with the molecule via its dipole moment and polarizabibility, thereby 
neglecting higher order contributions due to (higher) hyperpolarizabilities. Thus, in the framework of the Born- 
Oppenheimer approximation the rovibrational Hamiltonian reads 

H = Tr + ^-^^-^ +V{R) - FD{R) cose - — [a^{R) sin^ + a\\{R) cos^ 0] , (1) 

where R and 0, (p are the internuclear distance and the Euler angles, respectively, and we use the molecule fixed frame 
with the coordinate origin at the center of mass of the nuclei. Tr is the vibrational kinetic energy, fi,J(0, 0) is the 
orbital angular momentum, /i is the reduced mass of the nuclei, and V{R) is the field-free electronic potential energy 
curve (PEC). The electric field is taken oriented along the 2;-axis of the laboratory frame with strength F. The last 
three terms provide the interaction between the electric field and the molecule via its permanent electronic dipole 
moment function (EDMF) D{R), and its polarizability, with a||(i?) and a±{R) being the polarizability components 
parallel and perpendicular to the molecular axis, respectively. 

In the presence of the electric field, the dissociation threshold changes and it is given by the quadratic Stark shift 
of the free atoms, i.e., Edt{F) — —0.5 * F^{ai + 02), with ai i = 1,2 being the polarizabilities of the free atoms. 
In the presence of the electric field, only the azimuthal symmetry of the molecular wavefunction holds and therefore 
the magnetic quantum number M is retained. In this work we focus on levels with vanishing magnetic quantum 
number M = 0. For reasons of addressability, we will label the electrically-dressed states by means of their field-free 
vibrational and rotational quantum numbers (v, J) . 

Let us briefly investigate under which conditions the contribution of the molecular polarizabilities can be neglected 
in the Hamiltonian ([T]). For simplicity reasons, and without loss of generality, we use the effective rotor approach 
[37| . assuming that the rotational and vibrational energy scales differ significantly and can therefore be separated 
adiabatically, and that the field influence on the vibrational motion is very small and can consequently be treated by 
perturbation theory. Then, in the framework of this approximation the rovibrational Hamiltonian ([T]) is reduced to 

- S.J^ - F{D)l cos0-^ [{a^)l + {Aa)l cos^ 0] + E^^^o, (2) 

where = jjiiR^^)'^ is the field-free rotational constant of the state with quantum numbers v, J = and M = 0, 
iP1qq{R) and -E^QQ are the vibrational wave function and energy, respectively, and we encounter the expectation values 
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{R-^)l = (€ool^-'l€oo), {D)l = (V°ool^(^)l€oo), = (^°ool«^(^)l€oo), and (Aa)0 = (V°ool«l|(^) 

a^(i?)|'0l^oo)- Within this approach, at a certain field strength F the interaction due to the polarizabihty can be 

neglected in the effective rotor Hamiltonian ([2]) if 



2{D)l 



>> 1 and 



2{D)l 
F(Aa)u 



>> 1. 



To analyze the very weak field regime, we rescale the effective rotor Hamiltonian ^ with , and assume that the 

ratios -^-{a^)^, and ■^-{Aa)° are smaller than the rescaled field-free rotational kinetic energy J{J + 1). 

Then, for a certain state with quantum numbers v, J and M, time-independent perturbation theory provides the 
following second order correction to the field-free energy 



E, 



(2) _ 
u,J,M ~ 



A 



JM- 



(3) 



where the angular coefhcients (38| are given by 

Am 



2J(J+f)(2J+l)(2J + 3) 
f 

^6 



for J > 



(4) 



and 



CjM — 



.f2 



J2 - M2 



{j + iy-M 

(2J + 1)(2 J + 3) ^ (2 J + f )(2J - 1) ■ 



(5) 



This second order correction to the rotational energy depends on the molecular system and the symmetry of 
the considered state through the expectation values (D)^, {a±)^ and {Aa)1 . In the perturbative regime, the 



polarizabihty terms can be neglected if 



>> 1 and 



>> 1. We have Cjm > 



coefficient becomes increasingly smaller than Cjm for increasing values of J; for example for J 



Ajm\, and the 
15 and M = 



; X 10"'' and Ci5,o = 0.5005, and for J ^ 15 ^15,15 = -8.859 x 10"'' and Ci 



3.030 X 10"^ As a 



consequence we encounter the situation that for high rotational excitations of certain molecular systems the interaction 
due to the molecular polarizabihty could be the dominant one. We emphasize that the above considerations are valid 
only for weak fields and within the effective rotor approach. 



III. RESULTS 



In the present work, we have performed a full rovibrational study of the influence of an external static electric 
field on the highly excited rovibrational states of the LiCs molecule. The PEC, EMDF and polarizabihty components 
of the electronic ground state of LiCs are plotted as a function of the internuclcar distance in Figures [Ija), 
and (b), respectively. For the PEC, we use the experimental data of ref.|39j, which includes for the long-range 
behaviour the van der Waals terms, — J2n=6 8 10 ^n/ , and an exchange energy term, —AWe~^^, see ref. 39] for 
the values of these parameters. The EDMF and polarizabilities are taken from semi-empirical calculations performed 
by the group of Dulieu [H, H^]. The EDMF is negative and its minimum is shifted by 1.4 oq with respect to the 
equilibrium internuclcar distance Re = 6.94 ao of the PEC. For the electronic ground state of the polar alkali dimers 
the long-range behaviour of the EDMF is given by D-j/lC [4l|, this function has been fitted to the theoretical data 
for R > 18.15 ao with 1)7 = —5 x 10~^ a.u. Regarding the polarizabihty, both components smoothly change as R 
is enhanced and a±{R) > a\\{R) for any R value. They satisfy that limfl^oo Oi±{R) = limij^oo Oi\\i3) — ^Li + acs, 
with the polarizabilities of the Li and Cs atoms a^i = 164.2 a.u. and acs — 401 a.u., respectively, [i^, E^. Thus, for 
i? > 26 ao the theoretical data were extrapolated by means of exponentially decreasing functions to match the constant 
value aLi + CiCs- For computational reasons, ai_{R) and a\\{R) are extrapolated for i? < 5 and 4 ao, respectively. 
Since this study is focused on highly excited levels lying close to the dissociation threshold, we are aware of the fact 
that our results strongly depend on the assumptions made for the long-range behaviour of D(R), a±{R) and ck||(i?), 
and on the extrapolations performed at short-range for a±{R) and a||(i?). However, let us remark that the overall 
behaviour and physical phenomena presented here remain unaltered as these parameters are altered. 

For the lowest rotational excitations within each vibrational band of LiCs, we have investigated and compared the 
field interactions with the dipole moment and polarizabihty presented in the previous section. Within perturbation 
theory, the interaction due to the molecular polarizabihty becomes comparable to the one due to the dipole moment 
only for the last two vibrational bands. Assuming that the effective rotor conditions are satisfied, the interaction with 
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FIG. 1: (a): Electronic potential curve (solid) and electric dipole moment functions (dashed), and (b): parallel Q!||(-R) (solid) 
and perpendicular a±{R) (dashed) components of the polarizability of the electronic ground state of the LiCs molecule. 

the polarizability can be neglected for those levels with i/ < 47 and 48 < < 52 if the field strength is smaller than 
10"'^ and 2 x 10^^ a.u., respectively. Whereas, for the vibrational bands ly — 53 and 54 both interactions possess the 
same order of magnitude for the much weaker fields F Q x 10"^ and 8 x 10~^ a.u., respectively. Furthermore, the 
absolute values of the quadratic Stark shifts of the atomic energies are larger than the binding energies of the last 
bound state for F > 10"^ a.u., which also justifies that the interaction with the polarizability has to be included in 
the present study. 

Here, we consider the highest rotational excitations (M = 0) for the last four vibrational bands, 51 < < 54, of 
LiCs with binding energies smaller than 0.28 cm~^. We focus on the strong field regime F = 10~^ — 3.4 x 10"'' a.u., 
i.e., F = 5.14 — 1747.6 kV/cm, which includes the experimentally accessible range of strong static fields and possibly 
quasistatic fields. We remark that such strong fields are considered to induce the below-described peculiar behaviour 
of these states. Most of the overall probability is located at the outermost hump of these states, i.e., in regions where 
the EDMF possesses small values and the polarizabilities are close to au + acs- Thus, strong fields are needed in order 
to observe a significant field-effect on these levels. At these field strengths the corresponding rovibrational dynamics 
cannot be described by means of the effective or (due to avoided crossings) even the adiabatic rotor approximations 
[53, H|, and, of course not by perturbation theory Hence, the two-dimensional Schrodinger equation associated 
to the nuclear Hamiltonian ([1]) has to be solved numerically. We do this by employing a hybrid computational method 
which combines discrete and basis-set techniques applied to the radial and angular coordinates, respectively 13, ^5^. 

Since in the presence of the field the dissociation threshold is E]jt{F) = —0.5 * F^{a]^i + acs), we define the 
energetical shift with respect to this dissociation threshold as e^j = E^^j{F) — Eot{F)^ with E^,j{F) being the 
energy of the [v, J) state at field strength F. Figure Ufa) shows these Stark shifts e^^j satisfying Si,,.] > —0.28 cm"'^ 
in the above-provided range of field strengths. This spectral window includes the (54,0), (54,1), (53,4), (53,5), 
(53,6), (52,10) and (51,15) states, and for F > 3 x 10~^ a.u. also the (52,9) level. Since all these levels possess 
the same symmetry (M = 0) for 7^ 0, and since the field strength is the only parameter at hand to vary the 
rovibrational energies, the von Neumann- Wigner noncrossing rule [45'] holds and we encounter only avoided crossings 
of energetically adjacent states but no exact crossings in the field-dressed spectrum. In the vicinity of the avoided 
crossing the nuclear dynamics is dominated by a strong interaction and mixing of the involved rovibrational states. 
For reasons of simplicity (according to the Landau-Zener theory) we assume that the avoided crossings are traversed 
diabatically as F is increased. Thus, a certain state has the same character before and after the avoided crossing, e.g. 
taking for example the (54, 0) level we observe that it keeps its high-field seeking trend, see Figure HJa). 

Before studying these avoided crossings in more detail, let us analyze the general behaviour of the binding energies. 
For all e^j we observe a very weak dependence on F for F < 3 x 10~^ a.u. The larger the field-free rotational quantum 
number of a state, the stronger is the field strength needed to encounter a deviation of e^j from its field-free value. 
With a further enhancement of the strength, e^^j increases (decreases) for the high (low)-field seekers. The strong field 
dynamics is dominated by pendular states, whose binding energies increase as F is augmented. Their main feature 
is their orientation along the field axis: They represent coherent superpositions of field- free rotational levels [38]. In 
our spectral region, this regime is only reached for the (54,0), (54,1) and (53,4) states. Indeed, £54,0 monotonically 
decreases as F increase, whereas £54. 1 and £53. 4 initially increase and reach broad maxima, decreasing thereafter. In 
contrast, the binding energies of the (52,9), (52, 10), (53,6) and (53,5) decrease as F is increased. Due to its large 
field-free angular momentum, the (51, 15) level is the least affected by the field. Initially, £51.15 is reduced as F is 
enhanced, passes through a broad minimum and increases thereafter. The contribution of the molecular polarizability 
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FIG. 2: (a) Energy shifts with respect to the dissociation threshold e^j, and expectation values (b) (coaO), (c) (J'^) and (d) 
(R), as a function of the field strength for the states with field-free vibrational and rotational quantum numbers (54, 0) (solid), 
(54,1) (solid thick), (53,4) (dotted-dashed), (53,5) (long dashed thick), (53,6) (long dashed), (52,10) (dotted), (52,9) (short 
dashed), and (51, 15) (double dotted-dashed). Note that in the panel (c) only the results for the levels (54,0), (54, 1), (53,4), 
(53, 5) and (53, 6) are included. 

causes that this state is a high-field seeker in the weak field regime, whereas if only the interaction with the dipole 
moment is included it has a low-field seeking character. 

Regarding the avoided crossings, some of them are very narrow and can not be identified as such on the scale 
used in Figure [2ja), e.g. those among the pairs of levels (54,1) — (53,6), (54,1) — (51,15), (54,0) — (51,15), and 
(54, 0) — (52, 10). In contrast, other avoided crossings are very broad, and they are characterized by a strong couphng 
between the involved molecular states. The avoided crossing between the (54, 1) and (53, 5) levels takes places at strong 
fields, and the minimal energetical gap is AE = \e54^ i — £53,5! = 1.13 x 10^^ cm^^ for F « 2.23 x lO^** a.u. The (54, 0) 
state is involved in an avoided crossing with the (53,5) level characterized by AE = |e54,o ~ £53,5! — 2.82 x 10^'^ 
cm~-^ for F w 1.399 x 10"'' a.u., and another one with the (53,4) state with AE = |£54, — £53,4! = 1-08 x 10"^ 
cm"-' for F « 1.918 x lO^"* a.u. The (53,4) level experiences an avoided crossing with the (52,9) state, with minimal 
energetical separation AE — |£53,4 — £52,9! = 1-12 x lO^^^ cm^^ for F w 3.25 x 10"'* a.u. For other alkali dimers weaker 
electric field strengths might suffice to exhibit similar avoided crossings. These avoided crossings are wide enough to be 
experimentally observed in a similar way as it has been done for the weakly bound spectrum of Cs2 dimer in a magnetic 
field |46j]. In this system couphng strengths, larger than 1.67 x 10^^ cm^' were experimentally estimated; i.e., 
even the energetical separation of the (54, 1) and (53, 6) avoided crossing, AE — [£54^ 1 — £53^ el = 9.6 x 10~^ cm~^ for 
F w 1.1 X 10"'' a.u., could be measured. Moreover, as it has been done for the Cs2 molecule [s^ in a magnetic 
field, a suitable electric-field ramp could be used to transfer population from high to low rotational excitations in a 
controlled way, by either diabatically jumping or adiabatically following these electrically induced avoided crossings. 

An interesting physical phenomenon is observed in the evolution of the (53, 6) level in the spectrum. £53^ g increases 
as F increases, and after passing the avoided crossing with the (54, 1), £53,6 becomes positive for F > 1.6 x 10~''a.u. 
The Stark increase of the (53, 6) energy surpasses the reduction of the dissociation threshold, and this level is shifted to 
the continuum. Hence, if the LiCs is initially in the (53, 6) level, it will dissociate as the field strength is adiabatically 
tuned and enhanced above F > 1.6 x lO"'* a.u. Therefore a channel for molecular dissociation is opened as the electric 
field is modified. Of course, the inverse process is also possible, and the continuum state formed by two free atoms can 
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be brought into a bound state by lowering the electric field strength. Indeed, it has been proved that a static electric 
field could be used to manipulate the interaction between two atoms such that a virtual state could be transformed 
in a new bound state, i.e., the molecular system supports a new bound level [23l |. 

To illustrate the appearance of this phenomenon for a low-lying rotational excitation, we have performed a similar 
study for a designed molecule. We have taken the theoretical PEC of the electronic ground state of LiCs 
computed by the group of AUouche [13] with the van der Waals long-range potential, Ce/R^, but modifying the LiCs 
Cq coefficient to Cg = 2225 a.u. The (54, 1) level is shifted towards the dissociation threshold, having a field-free energy 
Ec,4i w —5.9 X 10~^ cm~^. As electric dipole moment function and polarizabilities we have used the corresponding 
functions of the LiCs molecule described above. The last most weakly bound states of this toy system have been 
studied in the presence of a static electric field, but for the sake of simplicity we discuss here only the results for the 
(54, 1) level. As F is enhanced £54. 1 increases, and becomes positive for > 5 x 10~^ a.u.; note that this field strength 
is much weaker than the above used ones. For a bound level, a further enhancement of the field would change its 
character, and its binding energy would increase as F is augmented. We have observed the same phenomenon for the 
(54, 1) state, which becomes bound again for F > 1.7 x 10^"* a.u., and 654. 1 decreases thereafter. The level has been 
captured by the nuclear potential demonstrating that the reverse process is possible. Starting with two free atoms 
with the correct internal symmetry, by adiabatically tuning the field the dimer is formed in a highly excited level. 

Due to negative sign of the EDMF, the main feature of the pendular regime (focusing again on LiCs) is the 
antiparallel orientation of the states along the field axis. The orientation can be estimated by the expectation value 
(cos 6'): The closer |(cos6')| is to one, the stronger is the orientation of the state along the field. Figure[2)Jb) illustrates 
the evolution of (cos 9) as the field strength is changing. The initial behaviour of (cos 6) for weak fields depends on the 
character of the corresponding level. For the (54, 0) state {cos 6) monotonically decreases as F is increased, it achieves 
the largest orientation with (cos 6') < —0.7 for F > 5 x 10^^ a.u., except in the proximity of avoided crossings. For 
the (54,1), (53,5) and (53,4) levels, {cos6) reach a broad maximum decreasing thereafter. The orientation of the 
(54, 1) and (53, 4) states becomes antiparallel for stronger fields. Not considering the proximity of an avoided crossing 
region, the (54, 1) state shows a significant orientation with {cos 9) < —0.4 for F > 1.31 x 10^"* a.u. The remainder 
of states keep a pinwheeling character, and (cos 9) increases as F is augmented. Since we have used the notation 
that the avoided crossings are traversed diabatically, a certain state {cos 9) reestablishes its increasing or decreasing 
trend once the avoided crossing has been passed. The smooth behaviour of {cos 9) is significantly distorted by the 
presence of these spectral features, where due to the strong mixing and interaction among the two involved states 
(cos 9) exhibits sharp and pronounced maxima and minima. For example, the avoided crossing among the (54, 0) and 
(53,5) levels, is characterized by the values (cos 0)54, = —0.235 and (cos 0)53, 5 = —0.152, for F — 1.399 x 10^^ a.u., 
compared to the results (008^)54,0 = —0.847 and (0030)53^5 = 0.436 obtained for F = 1.3 x 10~^ a.u. Note that for 
the (54, 0) level (cos 6*) shows an additional maximum for F >3x 10^*, i.e., this level suffers another avoided crossing 
which is not observed in Figure[2{a), because £54,0 < —0.28 cm~^ for F > 2.54 x 10~^ a.u. 

The expectation value (J^) of the states (54,0), (54,1), (53,4), (53,5) and (53,6), is presented as a function of 
the electric field in Figurc[2fc). To provide a reasonable scale, the results for the (52, 10), (52, 9) and (51, 15) levels 
have not been included. This quantity provides a measure for the mixture of field-free states with different rotational 
quantum numbers J but the same value for M , i.e., it describes the hybridization of the field-free rotational motion. 
Analogous to the binding energy, (J^) shows for weak fields a plateau-like behaviour: The hybridization of the angular 
motion is very small and the dynamics is dominated by the field-free rotational quantum number of the corresponding 
state. For stronger fields, these states possess a rich rotational dynamics, with significant contributions of different 
partial waves, and (J^) decreases (increases) for the low-(high)-field seekers as F is enhanced. In the strong field 
regime, (J^) shows a broad minimum for the (54,1), (53,4) and (53,5) states, increasing thereafter. The pendular 
limit is characterized by the augment of (J^) due to the contribution of higher field- free rotational states. This regime 
is only achieved by the (54,0), (54, 1) and (53,4) levels. In contrast, the mixing with lower rotational excitations is 
dominant for the (53, 5) and (53, 6) states, and (J^) < J( J -|- 1), with J being the corresponding field-free rotational 
quantum number; similar results are obtained for the (52,9), (52, 10) and (51, 15) states not included in Figure [SJc). 
The presence of the avoided crossings significantly distorts the smooth behaviour of (J^). The (J^) of the level in 
an avoided crossing with the lowest (highest) field-free J exhibits a pronounced and narrow maximum (minimum) 
on these irregular regions. At the smallest energetical gap, we encounter similar values of (J^) for both states. For 
example, for F = 1.399 x lO"'' a.u. we obtain (J^) = 10.38 fi,^ and 12.40 fi^ for the (54, 0) and (53, 5) levels, respectively, 
compared to the values (J^)54,o = 3.04^12 and (J^)53,5 = 20.75 ^i^ for = 1.3 x 10""* a.u. 

The expectation value of the radial coordinate {R) is presented for these states and range of field strengths in Figure 
Wi<i). Only if the vibrational motion is affected by the field {R) should differ from its field- free value. Analogously 
to ejM and (J^), {R) represents approximately a constant for weak fields, and strong fields are needed to observe 
significant deviations from its field-free value. Indeed, the larger is the rotational quantum number of a state for 
F — 0, the least affected by the field is its {R). For the (54,0) level, (i?) monotonically decreases from 50.24 oq 
to 28.13 oq as F is enhanced from to 3.4 x 10^^ a.u. For the (54,1), (53,5) and (53,4) states, (i?) increases as 
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FIG. 3: Probability densities (a) of the state (54, 0) and (b) of the state (53, 5) for F = 1.3 x 10 a.u., i.e., close to the avoided 
crossing but still without mixing of the rovibrational field-dressed states. 



F is augmented, reaches a broad maximum and decreases thereafter. The (54, 1) level is significantly affected with 
a reduction from {R) — 52.52 ao to 30.96 ao for F = and 3.4 x 10"'' a.u., respectively. For the (53,4) state this 
effect is much smaller, and (i?) is modified from the field-free result 30.78 aq to 28.32 oq. For (53,5) we observe that 
for F — 3.4 X 10"'' a.u. (i?) is by 4.57 ao larger than its value for F = 0. (i?) increases as F is enhanced for the 
remaining states, their total rise being smaller than 5ao for the analyzed levels j/ = 51 and 52. As the (53,6) state 
is shifted to the continuum, the slope of (i?) becomes very steep, and {R) is enhanced from (i?) = 32.85 ao up to 
41.30 ao for F = and 1.6 x 10~^ a.u., respectively. The field effect on the vibrational motion can be explained 
as follows: The probability density of those levels with a antiparallcl (parallel) orientation is mostly located in the 
7r/2 <0<7r(O<6'< 7r/2) region, where the dipole moment interaction is attractive (repulsive). As a consequence, 
the wavefunctions are squeezed (stretched) compared to their field-free counterparts to reduce the energy. Again, in 
the vicinity of the avoided crossing {R) exhibits very similar values for the two involved states. For example, we have 
found that (i?) = 33.02 and 33.13 ao for the states (54,0) and (53,5) and F = 1.399 x 10~^ a.u., respectively. 

To gain a deeper insight into the coupling of the vibrational and rotational motions induced by the electric field, 
we have analyzed the corresponding wavefunctions of two states involved in an avoided crossing. As an example, 
we discuss here the (54, 0)-(53, 5) avoided crossing. For comparison, let us first analyze their wavefunctions for 
F — 1.3 X 10~^ a.u., i.e., 'below' the avoided crossing where the mixing is not yet appreciable. The contour plots of 
the probability densities, \iIj{R,9)\^ sin 6, in the [R^O) plane are presented in Figures EJ^a) and (b) for the (54,0) and 
(53, 5) states, respectively. Since most of the overall probability of these weakly bound levels is located in the outer 
most hump, the radial coordinate has been restricted in these plots to the interval 15ao < i? < 52 oq. Indeed, more 
than 89% of the (54,0) and (53,5) probability densities are located for R > 30 ao, and 25 ao, respectively. Due to 
the pronounced antiparallel orientation of the (54, 0) level, (cos 9) = —0.847, the corresponding probability density 
shows a pendular-like structure, it is located in the region 37r/4 < 6* < tt and the maximal value is obtained at 
— 2.77 and R = 35.87 ao. The typical oscillator-like behaviour with 6 maxima reminiscent from its field free angular 
momentum J = 5 is observed in the (53, 5) probability density, see Figure[3I^b). Since this state has still a pinwheeling 
character, the corresponding probability density is distributed over the complete interval < 9 < tt, however, due to 
the parallel orientation of this state, {cos 9) = 0.437, the probability density is larger in the region 9 < ^. Moreover, 
the influence of the field on the vibrational motion provokes an inclination of the intcrnuclear axis of this level, i.e., 
the corresponding wavefunction is stretched and squeezed in the regions 9 < ^ and 6* > f , respectively. The squeezing 
effect also appears for the (54,0) state which possesses a strong antiparallel orientation (see Figure [Sja)). 

As the electric field is enhanced approaching the region of the avoided crossing, a strong interaction between the 
involved states takes places and the rovibrational dynamics is affected drastically. The contour plots of the (54, 0) 
and (53,5) states for F = 1.399 x 10"'* a.u., which corresponds with the minimal energetical gap between them, 
are shown in Figures [4{a) and (b), respectively. Although, at this field strength their orientation and hybridization 
of the angular motion are very similar, there exist significant differences with respect to their wavefunctions. The 
above-described regular structures typical for an oscillator and pendular-like distributions are lost. Even more, for 
both levels it is not possible to identify an orientation of the molecule, and the most pronounced maxima are not 
necessarily located at the outermost turning points. The (54, 0) probability density is distributed in the interval 
< 9 < TT, the largest maximum is at = 2.701 and R = 34.02 ao, and it is accompanied by several less pronounced 
maxima at smaller 9 values. The (53, 5) probability density exhibits three maxima with similar probability density, 
the first one at 6* = 0.44 and R = 35.24 ao, the second one a.t 9 = 2.95 and R — 33.73 ao, and the third one, at 
9 = 2.70 and R = 25.97 ao, which is shifted towards smaller internuclear separations from the outermost turning 
point. Both configurations exhibit significantly distorted patterns, and they show a strong mixing between the radial 
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FIG. 4: Probability densities (a) of the state (54,0) and (b) of the state (53,5) at the field strength F = 1.399 x 10"" a.u. 



and angular degrees of freedom. In general, at the avoided crossings the nuclear dynamics of the field-dressed states 
are characterized by an asymmetric and strongly distorted behavior, exhibiting pronounced localization phenomena. 

We have also analyzed these weakly bound levels taking into account only the interaction of the field with the 
permanent electric dipole, and not considering the contribution of the molecular polarizability. The results look 
qualitatively similar but show a quantitatively different behaviour as a function of the electric field, and the effect 
of the polarizability becomes important for F > 10~^ a.u. The polarizability terms cause the mixing of states with 
field-free rotational quantum numbers J and J ± 2. If the polarizability is included, the dissociation energies are 
smaller, i.e., for a certain F value the modulus of the displacement of the dissociation threshold is larger than the 
modulus of the energetical shift due to polarizability of a certain level, and the avoided crossings are also broader. 
Without the contribution of the polarizability term the (53, 6) level is not shifted to the continuum as the field is 
increased, and the inverse phenomenon appears, i.e., the (54, 2) level becomes a bound state for F > 1.8 • 10~* a.u. 

The validity of the adiabatic and effective rotor approaches has been previously demonstrated for vibrational low- 
lying levels of the LiCs dimer [33| . However, this does not hold true for the part of the spectrum considered here. Since 
both approximations do not include the full coupling between the vibrational and rotational motion, the presence 
of the avoided crossings in the spectrum is not reproduced. In addition, significant errors are found for the binding 
energies and the expectation value (i?) of the rotational excitations even in the absence of the field, e.g. the (53, 6) 
and (51, 15) levels are not bound within these approaches. The above-discussed field-effects on the vibrational motion 
can not be explained using an effective rotor description [s^]- However, the adiabatic results qualitatively reproduce 
the orientation and hybridization of the angular motion as well as the stretching and squeezing of the vibrational 
motion. Numerically significant deviations are encountered in the avoided crossing regions. 



IV. CONCLUSIONS AND OUTLOOK 



We have investigated the influence of a strong static and homogeneous electric field on the highly excited rovibra- 
tional states of the electronic ground state X^I]+ of the alkali dimer LiCs by solving the fully coupled rovibrational 
Schrodinger equation. The interaction of the field with the electric dipole moment function as well as with the 
molecular polarizability has been taken into account. We focus here on the last rotational excitations with vanishing 
azimuthal symmetry within the last four vibrational bands, 51 < < 54. Due to their large extension, strong fields 
are needed in order to observe a significant field influence. 

The richness and variety of the resulting field-dressed rotational dynamics has been illustrated by analyzing the 
energetical Stark shifts, as well as, the orientation, the hybridization of the angular motion and the vibrational 
stretching and squeezing effects. Whether we encounter a squeezing or stretching of the vibrational motion depends 
on the angular configuration: The molecule tries to minimize its energy leading to stretching for a parallel configuration 
and squeezing for an antiparallel one. In the strong field regime, the electrically-dressed spectrum is characterized by 
the presence of pronounced avoided crossings between energetically adjacent levels. These irregular features lead to 
a strong field-induced mixing and interaction between the states, and they cause strongly distorted and asymmetric 
features of the corresponding probability densities. We stress the importance of identifying these irregular features: 
Their presence affects the radiative decay properties of the dimer, such as lifetime and transition probability for 
spontaneous decay, and they might significantly alter the chemical reaction dynamics. Even more, one of their 
possible applications is their use to transfer population between the involved states. 

We have demonstrated that if the last most weakly bound state is a low-field seeker it is possible to shift it to the 
atomic continuum by tuning the electric field, i.e., the molecular system dissociates into free atoms. The reverse process 
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is also possible, i.e., a continuum state, formed by two free atoms with the correct field- free rotational symmetry, 
can be transferred to a weakly bound molecular state by changing the field strength. These results suggest that by 
properly increasing or decreasing the value of the electric field, one could study in a control way the opening of a 
dissociation or an association channel. 

Although our study is restricted to a LiCs dimer and to the spectral region close to the dissociation threshold of its 
electronic ground state, we stress that the above-observed physical phenomena are expected to occur in many other 
polar molecules. 
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